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ABSTRACT 

Researchers  and  scientists  have  been  studying  Aloha  networks  for  many  years  hoping 
to  improve  the  channel  throughput  and  overcome  the  inherent  instability.  In  this  study,  we 
are  going  to  describe  some  of  the  Aloha  network  features.  An  application  of  the  slotted 
Aloha  network  is  considered  for  local/mobile  radio  communications.  The  near/far  effect, 
Rayleigh  fading,  and  shadow  fading  arise  in  local/mobile  radio  communications  due  to 
differences  in  distances  and  topography  between  each  of  the  users  and  the  base  station. 
These  various  effects  cause  a  packet  capture  effect  which  improves  overall  channel 
throughput  but  leads  to  different  packet  delay  times  for  the  various  users.  An  analysis  of  a 
ring  model  for  local/mobile  radio  communications  with  correct  packet  capture  is  considered 
in  this  paper.  The  correct  packet  capture  effect  of  one-ring  and  three-ring  networks  are 
studied. 

A  Markov  model  is  developed  for  a  slotted  Aloha  network  with  capture.  It  is  shown  that 
the  throughput  in  such  a  network  is  markedly  greater  than  'he  standard  1/e.  Perhaps  even 
more  important  is  the  result  that  such  networks  are  more  stable  under  overload. 
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I.   INTRODUCTION 

An  analysis  of  the  performance  of  a  local/mobile  radio  communications 
system  utilizing  the  slotted  ALOHA  random  access  protocol  is  investigated  in  this 
thesis. 

A.       ALOHA  NETWORKS 

1.  Pure  Aloha 

In  Aloha  networks,  the  information  comes  in  packets.  In  a  conventional 
pure  ALOHA  network  the  users  transmit  their  packets  whenever  new  data  are  ready 
for  transmission.  When  two  or  more  users  transmit  their  packets  simultaneously,  a 
collision  will  occur  in  the  channel,  and  the  packets  may  overlap  partially  or 
completely  (Figure  1).  The  collision  destroys  the  information  in  the  packets,  and  the 
affected  users  have  to  retransmit  their  packets.  These  collisions  cause  a  major 
degradation  in  performance,  which  limits  the  maximum  channel  throughput  of  pure 
ALOHA  to  l/2e  or  0.184  [Ref.  1]. 

2.  Slotted  Aloha 

In  a  slotted  ALOHA  network,  time  is  segmented  into  slots.  Each  slot  has 
a  duration  equal  to  the  packet  transmission  time  or  length.  Users  are  required  to 
synchronize  the  transmission  of  their  packets  so  that  the  leading  edges  of  the 
packets  are  aligned  with  the  beginning  of  the  predetermined  time  slot  at  the  receiver 
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tagged  packet 

packet 

Figure  1. 


The  Collision  in  Pure  ALOHA  Network. 


Since  packets  arrival  at  the  receiver  are  synchronized,  when  more  than  one  packet 
is  transmitted  simultaneously,  a  collision  will  occur  in  the  channel,  and  the  packets 
will  overlap  completely  (Figure  2).  The  elimination  of  partial  packet  overlap 
substantially  improves  the  performance  of  the  network.  The  maximum  channel 
throughput  is  increased  to  1/e  or  0.368  (twice  the  maximum  channel  throughput  of 
pure  Aloha)  [Ref.  1].  This  maximum  channel  throughput  is  not  necessarily  the 
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Figure  2.  The  Collision  in  Slotted  ALOHA  Network, 

throughput  attained  by  the  system.  In  heavy  traffic,  the  likelihood  of  collisions 
increases  significantly.  After  each  collision,  the  affected  users  have  to  retransmit  their 
collided  packets.  As  the  number  of  retransmitted  packets  becomes  large,  the  channel 
will  reach  a  point  where  a  collision  occurs  every  time  slot.  Hence,  no  packets  can 
be  transmitted  successfully  and  the  system  is  said  to  be  unstable. 


B.  ALOHA  NETWORK  IN  LOCAL/MOBILE  COMMUNICATIONS 

The  Aloha  network  is  receiving  a  great  deal  of  interest  for  use  in  applications 
which  involve  large  physical  areas  (eg.  ships,  multi-story  buildings,  warehouses,  urban 
environments,  etc.)  and  require  random  multiple  access.  The  two  primary  limitations 
of  the  Aloha  network  are  its  poor  maximum  throughput  and  its  inherent  instability. 
However,  the  Aloha  protocol  is  still  desirable  because  it  is  easy  to  implement. 

C.  PACKET  CAPTURE  EFFECT 

1.  Advantages 

In  local/mobile  radio  communications,  each  transmitter  produces  a 
different  signal  power  level  at  the  receiver  due  to  the  differences  in  receiver- 
transmitter  range,  shadow  fading,  and  Rayleigh  fading.  The  arriving  packet  with  the 
highest  power  level  will  have  the  best  chance  to  capture  the  receiver.  As  a  result,  it 
is  possible  for  the  receiver  to  accurately  decode  one  of  the  packets  involved  in  a 
collision.  Since  the  collision  does  not  necessarily  destroy  all  involved  packets,  the 
channel  throughput  will  increase  considerably.  The  capture  effect  significantly 
improves  the  performance  of  the  Aloha  network  [Ref.  2]. 

2.  Limitations 

The  advantage  of  improving  channel  throughput  is  achieved  in  favor  of 
users  near  the  base  station  and  at  the  expense  of  users  far  from  the  base  station.  The 
user  far  from  the  base  station  will  have  a  low  probability  of  capturing  the  base 


station  receiver,  and  consequently  will  have  to  transmit  its  packet  more  often  than 
a  user  near  the  base  station  and  hence  suffer  a  long  packet  delay  [Ref.  2]. 

D.  PERFORMANCE  MEASUREMENT 

Several  methods  have  been  considered  to  measure  the  performance  of  the 
Aloha  system.  They  are  packet  reception  probability,  the  correct  reception  of  the 
packet  address  header,  and  capture  measures  based  on  a  threshold.  Since  the  packet 
reception  probability  provides  the  best  measure,  this  method  is  selected  to  measure 
the  performance  of  the  system.  A  transmitted  packet  is  considered  to  have  captured 
the  base  station  receiver  only  when  the  entire  packet  is  received  correctly  in  the 
presence  of  one  or  more  interfering  packets  [Ref.  3]. 

E.  ORGANIZATION  OF  THE  THESIS 

The  analysis  is  divided  into  four  sections.  Section  II  investigates  the  capture 
effect  for  a  ring  model  in  local/mobile  communications.  Section  III  presents  the 
numerical  results  for  a  three-ring  network.  Section  IV  analyzes  the  system  throughput 
for  a  multiple  ring  Aloha  network  using  a  Markov  chain  model.  Finally,  discussions 
and  conclusions  are  presented  in  Section  V. 


II.   PACKET  CAPTURE  PROBABILITY  OF  A  RING  MODEL 

The  analysis  of  the  packet  capture  effect  in  a  ring  model  for  local  /mobile 
radio  systems  will  be  presented  in  this  section.  The  probability  of  correct  packet 
capture  will  be  evaluated  for  a  system  using  ideal  coherent  binary  phase-shift-keying 
(BPSK)  modulation. 

A.       SYSTEM  MODEL 

A  network  consisting  of  uniformly  distributed  users  in  an  annular  ring  will 
be  studied.  A  uniformly  distributed  density  function  for  the  user  traffic  is 
assumed.  The  model  illustrated  in  Figure  3  is  the  system  model  for  this  study. 

The  traffic  generated  within  a  distance  r  of  the  base  station  is 


G(r)  =  2ufxg(x)  dx  (1) 


where  r  is  the  distance  from  the  user  to  the  base  stati   ;  The  user  traffic  density  is 

g(r). 

The  total  offered  traffic  in  the  network  is  obtained     as  G   =   G(r  -»  »). 
Therefore 


G  =  2nfxg(x)dx  (2) 


Figure  3.  System  Model 


The  spatial  distribution  function  for  packet  generation  (  the  probability  of  a 
a  packet  generated  within  a  distance  R  <  r)  is  [Ref.  3] 


2tc 


FR(r)  =  PriRzr)  =  —fxg(x)dx 


(3) 


Therefore  the  probability  density  function  for  packet  generation  is  found  to  be 


A  W- 


rg(r) 


(4) 


The  probability  density  function  for  packet  generation  in  ring  R;,  defined  as 


ri<r<ri+1,  is 


aw= 


AW 


0  ,  otherwise. 


,  rfsr*rw 


(5) 


Assuming  that  users  are  uniformly  distributed  in  an  annular  ring  of  inner  radius 
Rj  and  outer  radius  R2,  we  have  the  user  traffic  density 


8(r)  = 


,  R^  r  z  R1 

n(R}-Rf) 

0,  otherwise 


(6) 


Substituting  (6)  into  (4),  we  have  the  probability  density  function  for  packet 
generation  as  follows 


f*(T)  = 


2r 


-,  R^r^ 
0  ,  otherwise 


(7) 


The  distribution  of  packet  generation  of  a  ring  with  Rj  <  r  <  R2  is  shown  in 
Figure  4.  The  probability  of  a  packet  being  generated  within  a  ring  of  radius  R  <  rt 
is  given  by 
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) 

k 

g(r) 

—    —    — 

R 

i                     R 

> 

2 

Figure  4.  User  Traffic  Density. 


w  ■  #/* 


2*  r         G 


G\    tc(/?22-/?,2) 


-dx 


(8) 


which  can  be  evaluated  to  yield 


2 

W  -  -r—2  (9) 

The  probability  of  a  packet  generated  within  a  ring  of  radius  R  <  ri+1  is 

*>,..)  ■  -£4  <10> 

R2-R1 

Substituting  (10)  into  (5)  and  (7),  we  get  the  probability  density  function  for 
packet  generation  in  ring  Rt 

where  i  equals  0,  1,  2,...,n-l  and  R;  =  {t1  <  r  <  ri+1} 

B.       BPSK  BIT  ERROR  PROBABILITY 

We  assume  that  users  communicate  with  the  base  station  using  coherent  binary 
phase  shift  keying  (BPSK).  For  ideal  coherent  binary  phase  shift  keying  (BPSK) 
modulation  in  a  Rayleigh  fading  channel,  the  received  signal  of  a  particular  user  0 
in  the  presence  of  n-1  interfering  users  is  given  by  (assuming  the  background 
Gaussian  noise  is  much  smaller  than  the  interference) 


'o(')  =coso(t)  +  J>ieAJ.(0  <12> 

i-1 
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where  Sj(t)  are  BPSK  signals,  the  C;'s  are  Rayleigh  distributed  signal  amplitudes,  and 
the  6j's  are  the  uniformly  distributed  phases  of  the  interfering  signals  (we  set  80  = 
0  and  assume  that  6;,  i  =  l,  2,...,  n-1  are  relative  to  this  reference). 
At  the  sampling  time  the  decision  variable  of  BPSK  receiver  is 

«-i 

vo  =  coao  +  Ec»  ai  cos6«  <13> 

1-1 

where  aj  =  {±1}  and  are  independent  and  identical  distributed  random  variables 
(i.i.d.).  The  second  term  on  the  right  side  of  the  above  expression  is  the  sum  of 
independent  Gaussian  random  variables. 
The  received  signal  power  is 

S  =  E  {(c0  atf)  =  E  {c02}  (14) 

The  received  noise  power  is 

N  =  E  {(v0-  c0  atf)  (15) 

Substituting  equation  (13)  into  equation  (15),  we  have  the  received  noise  power 
as  follows 


n-l 

N  =  E  {(£c.  a.  cos0()2)  (16) 

i-l 


Noting  that  the  mean  of  cos28j  is  1/2,  we  find 
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n-1 


L    i-1 


(17) 


Hence,  the  received  signal  to  noise  ratio  is 


S         IE  fc02) 
SNR  =  -  = 


N 


i-l 


(18) 


The  bit  error  probability  of  user  0  in  the  presence  of  n-1  interfering  users  for 
BPSK  in  a  Rayleigh  fading  channel  is  [Ref.  4] 


*       2 


1 

_ 

1 

\ 

u  2 

SNR 

(19) 


The  near/far  effect  is  modeled  by  assuming  that  the  received  power  w  conforms 
to  an  inverse  a  power  propagation  law;  that  is 


w  =  cr 


(20) 


where  r  is  the  distance  of  the  user  from  the  base  station  and  a  is  the  propagation 
constant.  Generally,  3<a<5  for  ultra  high  frequen  „  UHF)  applications.  The 
constant  c  depends  on  the  gains,  heights,  and  power  of  the  transmitter  and  receiver 
antennas.  Without  loss  of  generality,  we  assume  that  c  is  normalized  to  unity  for  all 
users;  that  is,  users  are  taken  to  be  identical  and  use  omnidirectional  antennas  with 
radiation  pattern  maxima  in  the  horizontal  plane. 
From  (20),  the  received  power  of  user  i  is 
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w,  -  E  ict }  =  rt~ 


Substituting  (21)  into  (18),  we  have 


SNR  = 


2wn  2rr 


n-l  n-l 


S>i    E'«"    E 


n-l  /  r  \-o 
i 


i=l 


i=l 


i=l 


vro/ 


and  the  bit  error  probability  is 


*       2 


1  - 

1 

n-l  /  r  \-a 

i+r, 

i 

\ 

<=i 

lroJ 

which  can  be  written  as 


1 

' 

1 

■ 

Pi.  = 

1   - 

fc 

2 

n-l  /  r   \a 

i+r, 

0 

N 

i.i 

lrJ 

(21) 


(22) 


(23) 


(24) 


For  numerical  analysis  purposes,  we  select  the  propagation  constant  a  =  4. 
Hence, 


(25) 


1 

1 

pb  - 

2 

1 

"-1  frA4 

i+r, 

0 

\ 

i=l 

Ui 

The  bit  error  probability  for  the  case  of  M-level  diversity  is  [Ref.  3] 
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(26) 


where  the  value  of  [i.  is  given  by 


V-  = 


i»l 


KriJ 


(27) 


C.      CODING  FOR  PACKET  CAPTURE 

Since  all  packets  are  transmitted  with  equal  power  and  because  of  the  near/far 
effect  and  the  Rayleigh  fading  channel,  the  operating  signal-to-noise  ratio  can  be  very 
low.  To  overcome  this  problem,  we  propose  a  combination  of  diversity  and  forward 
error  correction.  The  Viterbi  decoding  scheme  is  employed.  This  scheme  proposes 
maximum  likelihood  decoding,  which  is  relatively  easy  to  implement  for  codes  with 
small  memory  orders.  A  binary  convolutional  code  with  rate  R  =  1/2  and  constraint 
length  7,  and  Viterbi  hard  decision  decoding  is  considered.  Define  the  input 
modulation  as  M-ary  and  the  output  modulation  as  Q-ary.  When  binary  coding  is 
used,  the  modulator  has  only  binary  input  (M  =  2).  Similarly,  when  binary 
demodulator  output  quantization  is  used  (Q  =  2),  the  decoder  has  only  binary 
output.  In  this  case,  the  demodulator  is  said  to  make  hard  decisions  (M  =  Q).  To 
simplify  the  implementation,  most  coded  digital  communication  systems  use  binary 
coding  with  hard-decision  decoding.  Since  errors  at  the  output  of  the  Viterbi  decoder 
occur  in  burst  even  if  the  errors  at  the  input  of  the  decoder  are  independent,  the 
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probability  of  correct  packet  reception  cannot  be  evaluated  by  assuming  independent 
errors.  It  has  been  shown  that  probability  of  correct  packet  reception  Pc  is  lower 
bounded  by  [Ref.  3] 


Pc  *  a-pF)L 


(28) 


where  L  is  the  packet  length  and  PE  is  the  union  bound  first  event  error  probability 
[Ref.  5].  For  hard  decision  decoding,  a  tight  upper  bound  on  PE  is  given  by 


Pe*I1*A 


(29) 


d=d. 


where  ad  is  the  number  of  code  words  of  Hamming  weight  d,  df  is  the  free  distance 
of  the  convolutional  code,  and  Pd  is  the  first  event  error  probability  of  a  code  word 
of  weight  d.  The  first  event  error  is  made  with  probability  [Ref.  5] 


P.  =  { 


d 

E 

€=(<*+ 1)/2 


'd)pLV-PbM)d-e 


,  d  odd 


i(i) '« (1"'«)* +  £  I  - 1  **•  a'PmY"  - d  even 


e  =  l*d/2 


(30) 


The  bound  of  (29)  can  be  further  simplified  by  noting  that  for  d  odd 


p<< 


d 

E 


/^\ 


\ej 


P?M   V~PbM)m 


bW 


(31) 


Expressions  involving  PbM  can  be  moved  out  of  the  summation  since  PbM  is 
independent  of  the  summation  index  e.  As  a  result,  the  upper  bound  for  Pd  is 
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Taking  the  lower  limit  of  e  at  0  instead  of  (d+ 1)/2,  we  obtain 


P*  <  P£  Q-Pyj"  Y, 


'  <d\ 


(33) 


Ve/ 


The  sum  can  now  be  evaluated  to  obtain  an  upper  bound  on  Pd  of 


Pd  <  2d  p£  a-P^f2  (34) 


d        *     '  bM   Vx    '  bki 

Similarly,  for  d  even  an  upper  bound  on  Pd  is  found  to  be 


p&  v-pjr +  E  h  p%  v-pbJm    <35> 

t-da*\\e) 


V  d 

Pd<- 


2 

This  expression  can  be  simplified  in  a  manner  similar  to  that  used  to  obtain 
(34)  to  get  an  upper  for  Pd  when  d  is  even 

Pd<2d  PSil-Pu/12  06) 

which  is  identical  to  (34). 

The   correct  packet  capture   probability   q(n)   is   obtain  by  averaging  the 
probability  of  correct  packet  reception  over  all  possible  location  of  n  users.  Thus 

q(n)  =  jj...jPcfR{r^Urn_x)dr0...drn_l  (37) 

where  fR(r)  is  given  by  (7).  Now  substituting  (30)  into  (42),  we  obtain  a  lower  bound 
for  the  correct  packet  capture  probability  as 
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*2      *2 


q(n)  ±  j...j{l-P^r{r^..fIfrn_l)dr0..uirn_l  (38) 


*i   *i 
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III.   NUMERICAL  RESULTS 

A  numerical  analysis  has  been  performed.  Programs  using  Matlab  provide 
numerical  results  for  the  above  correct  packet  capture  probability  for  user  0.  Eight 
interfering  packets  have  been  used  in  the  analysis.  As  the  number  of  the  interfering 
packets  increases  the  computation  time  becomes  very  lengthy.  It  took  approximately 
24  hours  to  compute  q(l),  q(2)...q(6).  Therefore,  the  simplified  Pdin  (36)  is  used  to 
calculate  q(l),  q(2),...  q(7),  q(8).  The  result  for  q(n)  using  the  simplified  bound  of  Pd 
in  (36)  yields  results  that  are  very  close  to  those  obtained  with  (30).  For  each 
additional  user  after  eight  interferers,  the  results  do  not  change  significantly.  The 
source  codes  in  Appendix  A  are  executed  to  obtained  the  results  which  are  listed  in 
Table  II,  III,  IV  ,  and  IV.  The  results  and  source  codes  using  (30)  are  listed  in 
Appendix  B  for  reference  purposes. 

Numerical  results  are  obtained  for  the  specific  case  of  L  =  400  bits  and  a  = 
4.  The  levels  of  diversity  considered  are  M  =  1,  2,  and  4.  A  convolutional  code  of 
rate  1/2,  constraint  length  7  is  used.  The  values  of  ad  and  d  are  listed  in  Table  I. 
[Ref.  3] 
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Table  I. 


THE    COEFFICIENTS 
OF  ad  AND  d 


ad 

d 

11 

10 

38 

12 

193 

14 

1,331 

16 

7,275 

18 

40,406 

20 

We  will  study  a  ring  model  for  the  following  cases. 

A.       A  ONE-RING  NETWORK 

A  one-ring  is  used  for  the  network.  The  user  r0  locates  in  a  range  from  R^  = 
0  km  to  R2  =  50  km,  any  other  users  will  vary  in  the  same  range  with  user  r0.  The 
probability  of  correct  packet  capture  for  one-ring  is 

5050     50 


<,(«)  =  jj..Ja-PEf/R(rMrl)...fl)(rn.l)dr0...drn_l       (39) 


0  0       0 


where 
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w- 


2r. 


±  ,  Ozr,  s50 
502  j 

0    ,  otherwise 


(40) 


where  j  =  0,  1,  2,  3,  ...,  n-1.  The  calculated  values  for  q(n)  are  listed  in  Table  II. 


Table  II.  SIMULATION    RESULTS    FOR    THE    PROBABILITY    OF 

CORRECT    PACKET    CAPTURE    q(n)    FOR     ONE-RING 
NETWORK 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

q(7) 

q(8) 

M  =  l 

1 

.2049 

.1116 

.0789 

.0587 

.0459 

.0382 

.0334 

M  =  2 

1 

.3676 

.2142 

.1508 

.1116 

.0876 

.0741 

.0649 

M  =  4 

1 

.6318 

.3964 

.2801 

.2132 

.1750 

.1464 

.1252 

B.       A  THREE-RING  NETWORK 

The  following  is  an  analysis  of  a  three-ring  model.  Ring  one  is  {0  km  <  r  <  25 
km};  ring  two  is  {25  km  <  r  <  38.5  km}  and  ring  three  is  {38.5  km  <  r  <  50  km}. 

1.      Ring  one 

The  probability  of  correct  packet  capture  for  ring  one  is 

50     50  25 


q(n)  =  /.../  f(l-P^L  fR(r0\RMrxyfR(^dro-drn^    (41> 


0       0     0 


where 


20 


2r„ 


,  Os  r0  *  25 


/AW  ::  1  252  '  "    ° 

0    ,  otherwise 


(42) 


2y 


AW 


,  0  <;  r;.  *  50 


502 
0  ,  otherwise 


(43) 


where  j  =  1,  2,  3  ...n-1.  The  calculated  values  of  q(n)  are  listed  in  Table  III. 

Table  III.  SIMULATION    RESULTS    FOR    THE    PROBABILITY    OF 

CORRECT  PACKET  CAPTURE  FOR  RING  ONE 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

q(7) 

q(8) 

M  =  l 

1 

.6758 

.4412 

.3221 

.2532 

.2080 

.1770 

.1543 

M=2 

1 

.8579 

.6960 

.5722 

.4801 

.4036 

.3396 

.3141 

M  =  4 

1 

.9042 

.8099 

.7197 

.6350 

.5580 

.4901 

.4615 

2.      Ring  two 

The  probability  of  correct  packet  capture  for  ring  two  is 

50     50  38.5 


q(n)  =/.../   j \l-Pj- Ur^AhYUrn-iWv^n-i    (44> 


0      0     25 


where 


WM 


2rn 


,  25  <;  r0  <;38.5 


38.52-252 

0       ,  otherwise 


(45) 
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w  -  < 


2r. 

— >-  ,  O*  r,  «j  50 
502  ' 


(46) 


0      ,otherwise 
where  j  =  1,  2,  3,...n-l.  The  calculated  values  of  q(n)  are  listed  in  Table  IV. 


Table  IV. 


SIMULATION    RESULTS    FOR    THE    PROBABILITY    OF 
CORRECT  PACKET  CAPTURE  FOR  RING  TWO 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

q(7) 

q(8) 

M  =  l 

1 

.0199 

0 

0 

0 

0 

0 

0 

M  =  2 

1 

.2897 

.0606 

.0100 

.0019 

0 

0 

0 

M  =  4 

1 

.4731 

.2581 

.1080 

.0446 

.0166 

.0063 

.005 

3.      Ring  Three 

The  probability  of  correct  packet  capture  for  ring  three  is 

50     50  50 


q{n)  =  /.../  {  d-PE)LfR(r0\R3VR(r1)...fR(rn_l)  (47) 


0      0  38.5 


where 


^ ,  38.5  <;  r    s  50 

:  \  502-38.52 

0  ,  otherwise 


(48) 
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w  - 


It. 

— -L  ,  0  <;  r,  <;  50 
502  ' 


(49) 


0    ,  otherwise 
where  j  =  0,  1,  2,...n-l.  The  calculated  values  of  q(n)  are  listed  in  Table  V. 


Table  V. 


SIMULATION    RESULTS    FOR    THE    PROBABILITY    OF 
CORRECT  PACKET  CAPTURE  FOR  RING  THREE 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

q(7) 

q(8) 

M  =  l 

1 

0 

0 

0 

0 

0 

0 

0 

M  =  2 

1 

.0133 

0 

0 

0 

0 

0 

0 

M  =  4 

1 

.2302 

.0210 

.0008 

0 

0 

0 

0 

Plots  of  q(n)  versus  n  for  case  1)  and  each  of  the  rings  for  case  2)  are  shown 
in  Figures  5,  6,  7,  and  8. 
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Probability  Of  Correct  Packet  Capture  For  One-Ring  Network 
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Probability  Of  Correct  Packet  Capture  For  Ring  One 
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Probability  Of  Correct  Packet  Capture  For  Ring  Two 
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Probability  Of  Correct  Packet  Capture  For  Ring  Three 
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IV.    SYSTEM  THROUGHPUT  FOR  A  THREE  RING  ALOHA  NETWORK 

The  analysis  of  a  ring  model  for  mobile/local  radio  communications  with 
correct  packet  capture  probability  has  been  performed  in  the  previous  section.  The 
numerical  results  show  that  users  who  are  located  close  to  the  base  station  are  more 
likely  to  capture  the  base  station.  The  throughput  for  each  ring  will  be  calculated  by 
using  a  Markov  chain  model.  The  Markov  model  provides  a  detailed  model  of  the 
network's  dynamic  behavior;  therefore,  it  is  a  desirable  model  for  the  analysis. 

The  probability  qj  that  one  of  i  colliding  packets  will  survive  a  collision  and  that 
the  other  i-1  packets  will  need  retransmission  later  has  been  calculated  in  Section  III. 
The  previous  system  model  of  an  Aloha  network  with  N  terminals  (users)  sharing  one 
channel  for  transmission  to  one  receiver  (base  station)  is  used.  It  is  assumed  that 
positive  acknowledgements  from  the  receiver  are  never  lost,  so  the  user  will  not 
retransmit  unnecessarily.  This  assumption  has  been  shown  to  be  very  realistic  since 
modern  mobile  networks  usually  have  to  guarantee  a  n,v  "ty  five  percent  coverage  of 
the  serving  area.  The  values  of  qj  are  calculated  in  Section  III,  and  it  is  assumed  that 
q;  does  not  equal  to  zero  for  i  greater  than  one.  This  assumption  also  depends  on  the 
specific  properties  of  the  network.  In  an  undisturbed  channel,  qj  is  equal  to  1  and  q0 
is  equal  to  0  [Ref.  8]. 
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During  any  time  slot  each  terminal  (user)  will  be  in  any  of  three  states:  the 
original  state  O,  the  transmission  state  T,  or  the  retransmission  state  R.  The  terminal 
will  be  in  state  O  when  there  is  no  packet  ready  to  transmit  or  if  a  new  packet  has 
just  arrived  while  it  is  waiting  for  the  next  time  slot.  The  terminal  will  be  in  state  T 
when  it  is  actually  transmitting  a  packet,  successfully  or  not.  The  terminal  will  be  in 
state  R  after  it  has  transmitted  a  packet  unsuccessfully  or  when  it  is  waiting  for  a 
future  time  slot  for  retransmission.  The  terminal  will  be  designated  as  an  O-,  T-,  or 
R-terminal  depending  on  which  state  it  is  in.  It  is  assumed  that  the  process  which 
generates  new  packets  will  be  halted  in  a  terminal  while  it  has  a  packet  to  transmit. 
The  possible  state  transitions  for  an  individual  terminal  are  shown  in  Figure  9. 

It  is  assumed  that  changes  of  state  can  occur  at  the  beginning  or  at  the  end  of 
the  time  slot  in  the  Aloha  network.  An  O-terminal  can  either  remain  in  state  O  (do 
nothing)  or  change  into  a  T-terminal  at  the  beginning  of  a  time  slot.  Upon  changing 
into  a  T-terminal,  it  starts  the  transmission  of  a  new  packet.  The  generation  rate 
which  a  terminal  changes  from  state  O  to  state  T  is  denoted  p0.  An  R-terminal  can 
either  stay  in  state  R  (do  nothing)  or  change  into  a  T-terminal  at  the  beginning  of 
the  time  slot.  In  the  latter  case,  it  starts  the  retransmission.  The  retransmission  rate 
when  a  terminal  changes  from  R  to  T  will  be  denoted  pr.  At  the  end  of  a  time  slot, 
a  T-terminal  can  either  go  to  O  after  a  successful  transmission  or  to  R  after  an 
unsuccessful  transmission. 
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Po  1-Qi 

do  nothing       (  O  T  R  )      do  nothing 


Figure  9.  The  Possible  State  Transition  For  Each  Individual  Terminal 

A  slotted  Aloha  network  can  be  described  as  a  homogenous  finite  Markov 
process  with  N+l  possible  states  (0,  1,  ...N)  corresponding  to  the  number  of  R- 
terminals  at  the  end  of  each  time  slot. 

Given  0  <  n  <  N  and  0  <  m  <  N,  let  nnm  denote  the  state  transition  probability 
that  the  network  will  move  from  state  n  to  m  in  one  step.  There  are  three  different 
types  of  transitions: 

Type  1,  in  which  the  state  decreases  more  than  one.  This  is  not  possible 
because  only  one  packet  will  be  transmitted  successfully  in  one  time  slot. 

Type  2,  in  which  the  state  decreases  by  one  with  probability  q.  This  will  happen 
if  and  only  if  there  is  no  O-terminal  and  i  R-terminals  change  into  state  T  for  1  <  i 
<n  (one  packet  is  successful). 
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Type  3,  the  state  moves  to  m  >  n.  This  will  happen  with  probability  qi+m.n+1  if 
and  only  if  m-n+  1  O-terminals  and  i  R-terminals  change  into  state  T  for  0  <  i  <  n; 
or  with  probability  l-qi+m.n  if  and  only  if  m-n  O-terminals  and  i  R-terminals  change 
into  state  T  for  0  <  i  <  n  (no  packet  successful). 

It  is  assumed  that  the  terminals  are  operating  independently.  The  generation 
rate  and  the  retransmission  rate  are  denoted  p0  and  pp  respectively.  The  probability 
of  correct  packet  capture,  qi?  was  obtained  in  Section  III.  The  value  of  nnm  is 
calculated  as  [Ref.  8] 


it. 


0 


n   /    \ 


i=i 


d-Po)^".  e     .  p;.  a-pjr1' *t 


\  */ 


n/n 


=    < 


(N-n 
N-m 


-a -p^ 


N-m-l         m-» 


m-n  \-> 

Po    L, 

i=Q 


/«\ 


\1/ 


.(l-p^-'.Pr 


f\       \  /i  \      N-m 

(1-/V-(1-W»)+- ~7'Po-  Qi+m-n+1 

m-n  +  i 


,  if  m  <  n-\ 
,  if  m  =  n-\ 

,  if  m  z  n 


(50) 


To  calculate  the  equilibrium  state  occupation  probabilities  irn\  an  arbitrary 
positive  number  is  assigned  to  nn*.  The  nn'  is  calculated  as  follows  [Ref.  6] 


*«  = 


71 


nji-l 


/  B-l 

It  m     1   ~   7         t  ,•  7C  j  , 

1-1     /  j      i      i,n-l 
1=0 


(51) 


Normalizing  the  wn\  one  gets 


it. 


it_  = 


N 


i=0 


(52) 
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The  expected  throughput  of  the  network  in  state  n  (the  probability  to  transmit 
a  packet  successfully  in  a  time  slot)  is  found  to  be 


N-n 


f.  -  E {"'" )  ■  pi  a-*r**.  e(")  •  p'r  •  a -*>"'• «.*         (53) 

The  throughput  for  a  ring  is 

«-£«„/.  (54) 

n=0 

The  system  throughput  is  shown  in  Figure  10.  The  curves  SI,  S2,  and  S3  are  the 
throughputs  of  ring  one,  two,  and  three,  respectively.  The  curve  S1  +  S2+S3  is  the 
overall  throughput  of  the  network.  The  source  code  for  Figure  10  is  listed  in 
Appendix  C. 
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System  Throughput  For  A  Three-Ring  Network 
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V.   CONCLUSIONS 

At  first,  an  introduction  to  Aloha  networks  have  been  made.  The  slotted  Aloha 
network  is  selected  in  this  study,  since  its  system  throughput  is  double  the  pure  Aloha 
system  throughput.  Then,  the  probability  of  correct  packet  capture  of  slotted  Aloha 
is  studied  for  mobile/local  radio  communications.  An  application  of  slotted  Aloha 
in  a  ring  model  has  been  developed.  In  this  model,  the  collisions  do  not  necessarily 
destroy  all  the  involved  packets.  With  the  help  of  this  model,  the  probability  of 
correct  packet  capture  has  been  investigated.  It  has  been  shown  that  for  a  ring 
network  the  system  throughput  is  greater  than  1/e,  and  the  ring  network  is  more 
stable  under  overload.  The  numerical  results  in  Section  III  have  been  plotted  in 
Figures  5,  6,  7,  and  8.  The  probability  of  correct  packet  capture  of  one-ring  and 
three-ring  networks  have  been  studied.  The  numerical  results  and  its  plot  in  Figure 
5  show  that  as  the  number  of  users  grows  from  one  to  six  the  probability  of  correct 
packet  capture  drops  rapidly.  For  each  additional  user  j  r  six  users,  the  probability 
of  correct  packet  capture  does  not  decrease  as  quickly  as  with  only  a  few  users. 
When  the  number  of  users  is  large,  the  system  reaches  a  steady  state.  The  probability 
of  correct  packet  capture  increases  as  M  (diversity  level)  increases. 

Numerical  results  for  a  three-ring  network  show  that  when  the  users  are  located 
in  the  inner  ring  (e.g.  ring  one)  they  will  have  a  better  chance  to  capture  the  receiver. 
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As  the  users  move  to  the  outer  ring  (eg.  ring  three),  the  probability  that  they  will 
capture  the  base  station  decreases  rapidly.  The  system  throughputs  for  each  ring  have 
been  calculated  in  Section  IV  based  on  the  probability  of  correct  packet  capture 
obtained  previously. 
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APPENDIX  A 

As  mentioned  in  Section  III,  equation  (30)  is  used  to  calculate  q(n)  in  Tables 
II,  III,  IV  ,  and  V.  The  source  codes  for  the  results  shown  in  Tables  II,  III,  IV,  and 
V  are  listed  below.  Since  q(6),  q(7),  and  q(8)  require  a  long  computation  time,  the 
program  is  written  to  calculate  each  value  of  q(n)  separately  for  large  value  of  n  (n 
>6). 

The  source  code  for  the  numerical  results  of  Table  II  is  listed  below.  The 
number  of  users  (nuser)  is  initialized  to  2,  3,...,  and  8  to  calculate  q(2),  q(3), ...,  q(8), 
respectively.  In  the  listed  program  the  maximum  number  of  users  (nuser)  is  equal  to 
eight.  For  a  one  ring  network,  Rj  and  R2  are  varied  from  0  km  to  50  km. 
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%(*  this  program  calculates  q(8)  for  one  ring  network  *) 

%(*  the  coefficient  values  for  d  and  ad  *) 

d  =[10  12  14  16  18  10]; 

ad  =  [11  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  Rx  and  R2  *) 
rl  =  [0  0]; 
r2  =  [50  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 

alpha =4; 

L=400; 

nuser=8; 

ml  =  [12  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l)  =  (r2(l)-rl(l))/interv; 

delta(2)  =  (r2(2)-rl(2))/interv; 

r(l,l)  =  rl(l)  +  delta(l)/2; 

r(2,l)=rl(2)  +  delta(2)/2; 

for  i=2:interv 

r(l,i) =r(l,i-l)  +  delta(l); 
end 
for  i=2:interv 

r(2,i)=r(2,i-l)  +  delta(2); 
end 

%(*  calculate  the  probability  density  function  fR  equation  (7)  *) 
for  i=l:interv 

fr(l,i)=2*r(l,i)/(r2(ir2-rl(ir2); 
end 
for  i  =  l:interv 

fr(2,i)=2*r(2,i)/(r2(ir2-rl(2)^2); 
end 
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%(*  this  loop  calculate  q(8)  for  all 
%value  of  M  =  l,  2,  and  4,  equation  (38)  *) 
for  12  =  1:3 
m=ml(i2); 
sum=0; 
for  a=l:interv 
for  b=l:interv 
for  c=l:interv 
for  f=l:interv 
for  g=l:interv 
for  h  =  l:interv 
for  k  =  l:interv 
for  o  =  l:interv 

index  =  [a  b  c  f  g  h  k  o]; 
muy=0; 


%(*  calculate  value  for  n,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy=muy+(r(l,index(l)/r(2,index(il  +  l)))Aalpha; 
end 
muy  =  1  /sqrt(  1  +  muy) ; 

%(*  calculate  value  for  Pb,  equation  (25)  *) 
pb  =  (l-muy)/2 

%(*  calculate  the  value  for  PbM,  equation  (26)  *) 

pbm  =  0; 

for  cnt  =  0:m-l 

pbm  =  pbm  +  (comb((m- 1  +  cnt),cnt)  *  ((( 1  +  muy)/2)Acnt)); 
end 
pbm  =  (pbAm)*pbm; 

%(*  calculate  the  value  for  Pd  for  d  even,  equqation  (36)  *) 
pd=  (2*sqrt(pbm*(l-pbm))).Ad; 

%(*  calculate  Pe,  equation  (29)  *) 

pe  =  ad*pd'; 

if  pe>l,pe  =  l;,end 
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%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 
inc  =  (l-pe)AL*fr(l,o)*fr(2,a)*fr(2,b)*fr(2,c)*fr(2,f);  ' 
inc  =  inc*fr(2,g)*fr(2,h)*fr(2,k); 
sum  =  sura +inc; 
end 
end 
end 
end 
end 
end 
end 
end 
sum  =  sum*delta(2)A(nuser-l)*delta(l); 
end 


The  source  code  for  the  numerical  results  of  Table  III  is  listed  below.  For  ring 
one,  Rx  varies  from  0  km  to  25  km,  and  R2  varies  from  0  km  to  50  km.  This  source 
code  is  used  to  calculate  q(8).  To  compute  q(2),  q(3),...,q(n),  the  number  of  users 
(nuser)  is  equal  to  2,  3,...,  n,  respectively. 
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%(*  this  program  calculates  q(8)  for  ring  one  *) 
%(*  the  coefficient  values  for  d  and  ad  *) 
d  =[10  12  14  16  18  10]; 
ad  =  [11  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  Ra  and  R2  *) 
rl  =  [0  0]; 
r2  =  [25  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 

alpha  =  4; 

L=400; 

nuser=8; 

ml  =  [12  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l)  =  (r2(l)-rl(l))/interv; 

delta(2)  =  (r2(2)-rl(2))/interv; 

r(l,l)  =  rl(l)+delta(l)/2; 

r(2,l)=rl(2)  +  delta(2)/2; 

for  i=2:interv 

r(l,i)=r(l,i-l)  +  delta(l); 
end 
for  i=2:interv 

r(2,i)  =  r(2,i-l)  +  delta(2); 
end 

%(*  calculate  the  probability  density  function  fR  eq  '.ion  (7)  *) 
for  i  =  l:interv 

fr(l,i)=2*r(l,i)/(r2(l)-2-rl(ir2); 
end 
for  i  =  l:interv 

fr(2,i)=2*r(2,i)/(r2(l)A2-rl(2)A2)j 
end 
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%(*  this  loop  calculates  q(8)  for  all 
%value  of  M  =  l,  2,  and  4,  equation  (38)  *) 
fori2  =  l:3 
m  =  ml(i2); 
sum=0; 
for  a=l:interv 
for  b=l:interv 
for  c=l:interv 
for  f=l:interv 
for  g=l:interv 
for  h  =  l:interv 
for  k  =  l:interv 
for  o=l:interv 

index  =  [a  b  c  f  g  h  k  o]; 
muy=0; 


%(*  calculate  value  for  /x,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy=muy+(r(l,index(l)/r(2,index(il  +  l)))Aalpha; 
end 
muy  =  l/sqrt(l  +  muy); 

%(*  calculate  value  for  Pb,  equation  (25)  *) 
pb  =  (l-muy)/2 

%(*  calculate  the  value  for  PbM,  equation  (26)  *) 

pbm  =  0; 

for  cnt=0:m-l 

pbm  =  pbm  +  (comb((m-l  +  cnt),cnt)  *  ((( 1  +  muy)/2)/vcnt)); 
end 
pbm  =  (pbAm)  *pbm; 

%(*  calculate  the  value  for  Pd  for  d  even,  equqation  (36)  *) 
pd  =  (2*sqrt(pbm*(l-pbm)))  "d; 

%(*  calculate  Pc,  equation  (29)  *) 

pe  =  ad*pd'; 

if  pe>l,pe  =  l;,end 
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%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 

inc  =  (l-pe)AL*fr(l,o)*fr(2,a)*fr(2,b)*fr(2,c)*fr(2,f); 
inc=  inc*fr(2,g)*fr(2,h)*fr(2,k); 
sum  =  sum+inc; 
end 
end 
end 
end 
end 
end 
end 
end 
sum  =  sum*delta(2)A(nuser-l)*delta(l); 
end 


The  source  code  for  the  numerical  results  of  Table  IV  is  listed  below.  For  ring 
two,  R2  varies  from  25  km  to  38.5  km,  and  R2  varies  from  0  km  to  50  km.  This 
program  is  used  to  calculate  q(8).  To  calculate  q(2),  q(3),...,  q(n),  the  number  of 
users  (nuser)  is  equal  to  2,  3,...,  n,  respectively. 
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%(*  this  program  will  calculate  q(8)  for  ring  two  •) 

%(*  the  coefficient  values  for  d  and  ad  *) 

d  =[10  12  14  16  18  10]; 

ad  =  [11  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  Rj  and  R2  *) 
rl  =  [25  0]; 
r2  =  [38.5  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 

alpha =4; 

L=400; 

nuser  =  8; 

ml  =  [12  4]; 

%(*  calculate  the  grid  points  of  integration  *) 
delta(l)  =  (r2(l)-rl(l))/interv; 
delta(2)  =  (r2(2)-rl(2))/interv; 
r(l,l)=rl(l)  +  delta(l)/2; 
r(2,l)=rl(2)  +  delta(2)/2; 
for  i=2:interv 

r(l,i) =r(l,i-l)  +  delta(l); 
end 
for  i=2:interv 

r(2,i)  =  r(2,i-l)  +  delta(2); 
end 

%(*  calculate  the  probability  density  function  fR  equation  (7)  *) 
for  i  =  l:interv 

fr(l,i)=2*r(l,i)/(r2(ir2-rl(ir2); 
end 
for  i  =  l:interv 

fr(2,i)=2*r(2,i)/(r2(l)^2-rl(2)-2); 
end 
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%(*  this  loop  calculates  q(8)  for  all 
%value  of  M  =  l,  2,  and  4,  equation  (38)  *) 
for  12=  1:3 
m=ml(i2); 
sum=0; 
for  a=l:interv 
for  b=l:interv 
for  c=l:interv 
for  f  =  l:interv 
for  g=l:interv 
for  h  =  l:interv 
for  k  =  l:interv 
for  o  =  l:interv 

index  =  [a  b  c  f  g  h  k  o]; 
muy=0; 


%(*  calculate  value  for  /x,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy = muy  +  (r(  l,index(  1  )/r(2,index(i  1  +  1 )))  Aalpha; 
end 
muy=l/sqrt(l  +  muy); 

%(*  calculate  value  for  Pb,  equation  (25)  *) 
pb  =  (l-muy)/2 

%(*  calculate  the  value  for  PbM,  equation  (26)  *) 

pbm  =  0; 

for  cnt  =  0:m-l 

pbm = pbm + (comb((m- 1  +  cnt),cnt)  *  ((( 1  +  muy  )/2)~cnt)) ; 
end 
pbm  =  (pb"m)  *  pbm; 

%(*  calculate  the  value  for  Pd  for  d  even,  equqation  (36)  *) 
pd  =  (2*sqrt(pbm*(l-pbm))).Ad; 

%(*  calculate  Pc,  equation  (29)  *) 

pe  =  ad*pd'; 

if  pe>l,pe  =  l;,end 
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%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 
inc  =  (l-pe)AL*fr(l,o)*fr(2,a)*£r(2,b)*fr(2,c)*fr(2,0; 
inc=inc*fr(2,g)*fr(2,h)*fr(2,k); 
sum=  sum  +  inc; 


end 
end 
end 
end 
end 
end 
end 
end 
sum  =  sum*delta(2)/v(nuser-l)*delta(l); 
end 


The  source  code  for  the  numerical  results  of  Table  V  is  listed  below.  For  ring 
three,  R2  varies  from  38.5  km  to  50  km,  and  R2  varies  from  0  km  to  50  km.  This 
program  is  used  to  calculate  q(8).  To  calculate  q(2),  q(3),...,  q(n),  the  number  of 
users  (nuser)  is  equal  to  2,  3,...,  n,  respectively. 
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%(*  this  program  will  calculate  q(8)  for  ring  three  *) 

%(*  the  coefficient  values  for  d  and  ad  *) 

d  =[10  12  14  16  18  10]; 

ad  =  [11  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  Rx  and  R2  *) 
rl  =  [38.5  0]; 
r2  =  [50  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 

alpha =4; 

L=400; 

nuser=8; 

ml  =  [12  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l)  =  (r2(l)-rl(l))/interv; 

delta(2)  =  (r2(2)-rl(2))/interv; 

r(l,l)  =  rl(l)  +  delta(l)/2; 

r(2,l)  =  rl(2)  +  delta(2)/2; 

for  i=2:interv 

r(l,i)=r(l,i-l)  +  delta(l); 
end 
for  i=2:interv 

r(2,i)=r(2,i-l)  +  delta(2); 
end 

%(*  calculate  the  probability  density  function  fR  eq  r  ion  (7)  *) 
for  i  =  l:interv 

fr(l,i)  =  2*r(l,i)/(r2(ir2-rl(ir2); 
end 
for  i  =  l:interv 

fr(2,i)=2*r(2,i)/(r2(l)^2-rl(2)-2); 
end 
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%(*  this  loop  calculates  q(8)  for  all 
%value  of  M  =  l,  2,  and  4,  equation  (38)  *) 
for  i2  =  l:3 
m  =  ml(i2); 
sum  =  0; 
for  a=l:interv 
for  b  =  l:interv 
for  c=l:interv 
for  f  =  l:interv 
for  g=l:interv 
for  h  =  l:interv 
for  k  =  l:interv 
for  o  =  l:interv 

index  =  [a  b  c  f  g  h  k  o]; 
muy=0; 


%(*  calculate  value  for  i±,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy  =  muy+(r(l,index(l)/r(2,index(il  +  l)))/valpha; 
end 
muy=  l/sqrt(l  +  muy); 

%(*  calculate  value  for  Pb,  equation  (25)  *) 
pb  =  (l-muy)/2 

%(*  calculate  the  value  for  PbM,  equation  (26)  *) 

pbm  =  0; 

for  cnt=0:m-l 

pbm = pbm + (comb((m-l  +  cnt),cnt)  *  ((( 1  +  muy  )/2)/vcnt)) ; 
end 
pbm  =  (pb^m)  *  pbm; 

%(*  calculate  the  value  for  Pd  for  d  even,  equqation  (36)  *) 
pd  =  (2*sqrt(pbm*(l-pbm)))Ad; 

%(*  calculate  Pe,  equation  (29)  *) 

pe  =  ad*pd'; 

if  pe>l,pe  =  l;,end 
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%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 
inc=(l-pe)AL*fr(l,o)*fr(2,a)*fr(2,b)*fr(2,c)*fr(2,f); 
inc=inc*fr(2,g)*fr(2,h)*fr(2,k); 
sum  =  sum+inc; 
end 
end 
end 
end 
end 
end 
end 
end 
sum=sum*delta(2)^(nuser-l)*delta(l); 
end 
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APPENDIX  B 

Equation  (36)  is  used  to  calculate  q(n)  since  computing  q(n)  for  n  >  6  from 
equation  (30)  is  prohibitively  time-consuming.  The  following  programs  use  (30)  for 
Pd  to  calculate  q(n).  The  following  programs  are  listed  for  reference  purposes.  The 
source  codes  and  the  results  are  listed  below.  Since  q(n)  for  n  >  6  takes  a  long  time 
to  calculate,  the  program  is  written  to  calculate  only  up  to  q(6). 

The  number  of  users  (nuser)  is  initialized  to  2,  3,...,  6  to  calculate  q(2),  q(3), 
...,  q(6),  respectively.  In  the  listed  program  the  number  of  users  (nuser)  is  equal  to 
six.  For  a  one-ring  network,  Rx  and  R2  vary  from  0km  to  50km.  The  source  code  and 
its  simulation  results  are  listed  below  and  in  Table  VI,  respectively. 
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%(*  this  program  calculates  q(6)  for  one  ring  network  *) 

%(*  the  coefficient  values  for  d  and  ad  *) 

d  =[10  12  14  16  18  10]; 

ad  =  [ll  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  Rt  and  R2  *) 
rl  =  [0  0]; 
r2  =  [50  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 

alpha =4; 

L=400; 

miser  =  8; 

ml  =  [l  2  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l)  =  (r2(l)-rl(l))/interv; 

delta(2)  =  (r2(2)-r  l(2))/interv; 

r(l,l)  =  rl(l)  +  delta(l)/2; 

r(2,l)=rl(2)  +  delta(2)/2; 

for  i=2:interv 

r(l,i)=r(l,i-l)  +  delta(l); 
end 
for  i=2:interv 

r(2,i)=r(2,i-l)  +  delta(2); 
end 

%(*  calculate  the  probability  density  function  fR  equation  (7)  *) 
for  i  =  l:interv 

fr(l,i)  =  2*r(l,i)/(r2(ir2-rl(ir2); 
end 
for  i  =  l:interv 

fr(2,i)=2«r(2,i)/(r2(l)-2-rl(2)-2); 
end 
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%(*  this  loop  calculate  q(6)  for  all 
%value  of  M  =  l,  2,  and  4,  equation  (38)  *) 
fori2=l:3 
m=ml(i2); 
sum  =  0; 
for  a=l:interv 
for  b  =  l:interv 
for  c  =  l:interv 
for  f=l:interv 
for  g=l:interv 
for  o  =  l:interv 

index  =  [a  b  c  f  g  o]; 
muy=0; 


%(*  calculate  value  for  /z,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy=muy+(r(l,index(l)/r(2,index(il  +  l)))^alpha; 
end 
muy  =  l/sqrt(  1  +  muy); 

%(*  calculate  value  for  Pb,  equation  (25)  *) 
pb  =  (l-muy)/2 

%(*  calculate  the  value  for  PbM,  equation  (26)  *) 

pbm  =  0; 

for  cnt  =  0:m-l 

pbm =pbm  +  (comb((m-l  +  cnt),cnt)*(((l  +  muy)/2)~cnt)); 
end 
pbm  =  (pb"m)*pbm; 

%(*  calculate  the  value  for  Pd  for  d  even,  equqation  (30)  *) 
for  q=l:6 

pd(q)  =  0.5*comb(d(q),d(q)/2)*pbm^(d(q)/2); 

pd(q)=pd(q)*(l-pbm)A(d(q)/2); 

fore  =  l  +  d(q)/2:d(q) 

pd(q)=pd(q)  +  (comb(d(q),e)*(pbm"e)*(l-pbm)"(d(q)-e)); 

end 
end 
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%(*  calculate  Pe,  equation  (29)  *) 
pe  =  ad*pd'; 
if  pe>l,  pe  =  l;,  end 

%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 

inc=(l-pe)"L*fr(l,o)*fr(2,a)*fr(2,b)*fr(2,c)*fr(2,f);  ' 
inc  =  inc*fr(2,g); 
sum=sum+inc; 
end 
end 
end 
end 
end 
end 
sum=sum*delta(2)^(nuser-l)*delta(l); 
end 


Table  VI. 


SIMULATION    RESULTS    FOR    THE    PROBABILITY    OF 
CORRECT  PACKET  CAPTURE  FOR  ONE-RING  NETWORK 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

M  =  l 

1 

.228 

.138 

.096 

.075 

.059 

M  =  2 

1 

.412 

.238 

.159 

.123 

.096 

M=4 

1 

.636 

.442 

.314 

.241 

.197 

The  source  code  and  the  simulation  results  are  listed  below  and  in  Table  VII, 
respectively.  For  ring  one  Rx  is  varied  from  0  km  to  25  km,  and  R2  is  varied  from  0 
km  to  50  km.  This  source  code  is  used  to  calculate  q(6).  To  compute  q(2), 
q(3),...,q(n),  the  number  of  users  (nuser)  is  2,3,. ..,n  respectively. 
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%(*  this  program  calculates  q(6)  for  ring  one  *) 
%(*  the  coefficient  values  for  d  and  ad  *) 
d  =[10  12  14  16  18  10]; 
ad  =  [11  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  Rt  and  R2  *) 
rl  =  [0  0]; 
r2  =  [25  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 

alpha =4; 

L=400; 

nuser=8; 

ml  =  [12  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l)  =  (r2(l)-rl(l))/interv; 

delta(2)  =  (r2(2)-rl(2))/interv; 

r(l,l)  =  rl(l)  +  delta(l)/2; 

r(2,l)=rl(2)  +  delta(2)/2; 

for  i=2:interv 

r(l,i)=r(l,i-l)  +  delta(l); 
end 
for  i=2:interv 

r(2,i)=r(2,i-l)  +  delta(2); 
end 

%(*  calculate  the  probability  density  function  fR  equation  (7)  *) 
for  i  =  l:interv 

fr(l,i)=2*r(l,i)/(r2(l)^2-rl(l)A2); 
end 
for  i  =  l:interv 

fr(2,i)=2*r(2,i)/(r2(l)"2-rl(2)A2); 
end 
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%(*  this  loop  calculates  q(6)  for  all 
%value  of  M  =  l,  2,  and  4,  equation  (38)  *) 
for  12=1:3 
m=ml(i2); 
sum=0; 
for  a=l:interv 
for  b  =  l:interv 
for  c=l:interv 
for  f =  lrinterv 
for  g=  lrinterv 
for  o  =  lrinterv 

index  =  [a  b  c  f  g  o]; 
muy=0; 


%(*  calculate  value  for  /j,  equation  (27)  *) 
for  il  =  lmuser-l 

muy=muy+(r(l,index(l)/r(2,index(il  +  l)))^alpha; 
end 
muy  =  l/sqrt(  1  +  muy); 

%(*  calculate  value  for  Pb,  equation  (25)  *) 
pb  =  (l-muy)/2 

%(*  calculate  the  value  for  PbM,  equation  (26)  *) 

pbm  =  0; 

for  cnt  =  0:m-l 

pbm = pbm  +  (comb((m- 1  +  cnt),cnt)  *  ((( 1  +  muy)/2)^cnt)); 
end 
pbm  =  (pb^m)  *pbm; 

%(*  calculate  the  value  for  Pd  for  d  even,  equqation  (30)  *) 
for  q=  1:6 

pd(q)  =  0.5*comb(d(q),d(q)/2)*pbm^(d(q)/2); 

pd(q)  =  pd(q)*(l-pbm)A(d(q)/2); 

fore  =  l+d(q)/2:d(q) 

pd(q)=pd(q)  +  (comb(d(q),e)*(pbm"e)*(l-pbm)A(d(q)-e)); 

end 
end 
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%(*  calculate  Pe,  equation  (29)  *) 
pe  =  ad*pd'; 
if  pe>  1,  pe  =  l;,  end 

%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 
inc=(l-pe)"L*fr(l,o)tfr(2,a)*fr(2,b)*fr(2,c)*fr(2,0;  ' 
inc=inc*fr(2,g); 
sum  =  sum  +  inc; 
end 
end 
end 
end 
end 
end 
sum=sum*delta(2)A(nuser-l)*delta(l); 
end 


Table  VII.  SIMULATION    RESULTS    FOR    THE    PROBABILITY    OF 

CORRECT  PACKET  CAPTURE  FOR  RING  ONE 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

M  =  l 

1 

.751 

.505 

.355 

.2838 

.238 

M  =  2 

1 

.881 

.752 

.624 

.5169 

.437 

M=4 

1 

.916 

.836 

.761 

.688 

.617 

The  source  code  and  the  simulation  results  are  listed  below  and  in  Table  VIII, 
respectively.  For  ring  two,  Ra  is  varied  from  25  km  to  38.5  km,  and  R2  is  varied  from 
0  km  to  50  km.  This  program  is  used  to  calculate  q(6).  To  calculate  q(2),  q(3),..., 
q(n),  the  number  of  users  (nuser)  is  2,  3,...,  n  respectively. 
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%(*  this  program  will  calculate  q(6)  for  ring  two  *) 

%(*  the  coefficient  values  for  d  and  ad  *) 

d  =[10  12  14  16  18  10]; 

ad  =  [11  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  Rx  and  R2  *) 
rl  =  [25  0]; 
r2=[38.5  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 

alpha =4; 

L=400; 

miser =8; 

ml  =  [12  4J; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l)  =  (r2(l)-rl(l))/interv; 

delta(2)  =  (r2(2)-rl(2))/interv; 

r(l,l)  =  rl(l)  +  delta(l)/2; 

r(2,l)=rl(2)  +  delta(2)/2; 

for  i=2:interv 

r(l,i)=r(l,i-l)  +  delta(l); 
end 
for  i=2:interv 

r(2,i)=r(2,i-l)  +  delta(2); 
end 

%(*  calculate  the  probability  density  function  fR  equation  (7)  *) 
for  i  =  l:interv 

fr(l,i)=2*r(l,i)/(r2(ir2-rl(ir2); 
end 
for  i  =  l:interv 

fr(2,i)=2*r(2,i)/(r2(l)^2-rl(2)^2); 
end 
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%(*  this  loop  calculates  q(6)  for  all 
%value  of  M  =  l,  2,  and  4,  equation  (38)  *) 
for  12  =  1:3 
m  =  ml(i2); 
sum=0; 
for  a=l:interv 
for  b=l:interv 
for  c=l:interv 
for  f=l:interv 
for  g=l:interv 
for  o  =  l:interv 

index  =  [a  b  c  f  g  o]; 
muy=0; 


%(*  calculate  value  for  /x,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy=muy+(r(l,index(l)/r(2,index(il  +  l)))^alpha; 
end 
muy  =  l/sqrt(l  +  muy); 

%(*  calculate  value  for  Pb,  equation  (25)  *) 
pb  =  (l-muy)/2 

%(*  calculate  the  value  for  PbM,  equation  (26)  *) 

pbm  =  0; 

for  cnt  =  0:m-l 

pbm = pbm  +  (comb((m-l  +  cnt),cnt)*(((l  +  muy)/2)^cnt)); 
end 
pbm  =  (pb"m)*pbm; 

%(*  calculate  the  value  for  Pd  for  d  even,  equqation  (30)  *) 
for  q  =  l:6 

pd(q)  =  0.5*comb(d(q),d(q)/2)*pbm"(d(q)/2); 

pd(q)=pd(q)*(l-pbmr(d(q)/2); 

fore  =  l  +  d(q)/2:d(q) 

pd(q)=pd(q)  +  (comb(d(q),e)*(pbm"e)*(l-pbm)^(d(q)-e)); 

end 
end 
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%(*  calculate  Pe,  equation  (29)  *) 

pe=ad*pd'; 

if  pe>l,  pe  =  l;,  end 

%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 
inc=(l-pe)/sL*fr(l,o)*fr(2,a)*fr(2Jb)*fr(2,c)*fr(2,f); 
inc  =  inc*fr(2,g); 
sum=sun+inc; 
end 
end 
end 
end 
end 
end 
sum  =  sum*delta(2)^(nuser-l)*delta(l); 


end 


Table  VIII.  SIMULATION    RESULTS    FOR    THE    PROBABILITY    OF 

CORRECT  PACKET  CAPTURE  FOR  RING  TWO 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

M  =  l 

1 

.055 

.0018 

0 

0 

0 

M=2 

1 

.33 

.1021 

.0227 

.0052 

.0011 

M  =  4 

1 

.487 

.2971 

.1534 

.065 

.0293 

The  source  code  and  the  simulation  results  are  listed  below  and  in  Table  IX, 
respectively.  For  ring  three,  Rj  is  varied  from  38.5  km  to  50  km,  and  R2  is  varied 
from  0  km  to  50  km.  This  program  is  used  to  calculate  q(6).  To  calculate  q(2), 
q(3),...,  q(n),  the  number  of  users  (nuser)  is  2,  3,...,  n  respectively. 
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%(*  this  program  will  calculate  q(6)  for  ring  three  *) 

%(*  the  coefficient  values  for  d  and  ad  *) 

d  =[10  12  14  16  18  10]; 

ad  =  [11  38  193  1331  7275  40406]; 

%(*  the  number  of  interval  for  integration  *) 
interv=5; 

%(*  the  values  for  Rx  and  R2  *) 
rl  =  [38.5  0]; 
r2  =  [50  50]; 

%(*  set  value  of  a,  L,  number  of  user,  and  ml  for  M  *) 

alpha =4; 

L=400; 

nuser=8; 

ml  =  [12  4]; 

%(*  calculate  the  grid  points  of  integration  *) 

delta(l)  =  (r2(l)-rl(l))/interv; 

delta(2)  =  (r2(2)-r  l(2))/interv; 

r(l,l)  =  rl(l)  +  delta(l)/2; 

r(2,l)=rl(2)  +  delta(2)/2; 

for  i=2:interv 

r(l,i)=r(l,i-l)  +  delta(l); 
end 
for  i=2:interv 

r(2,i)=r(2,i-l)  +  delta(2); 
end 

%(*  calculate  the  probability  density  function  fR  equation  (7)  *) 
for  i  =  l:interv 

fr(l,i)=2*r(l,i)/(r2(ir2-rl(l)A2); 
end 
for  i  =  l:interv 

fr(2,i)=2*r(2,i)/(r2(l)"2-rl(2)"2); 
end 
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%(*  this  loop  calculates  q(6)  for  all 
%value  of  M  =  l,  2,  and  4,  equation  (38)  *) 
for  12  =  1:3 
m=ml(i2); 
sum=0; 
for  a=l:interv 
for  b=l:interv 
for  c=  l:interv 
for  f  =  liinterv 
for  g= liinterv 
for  o  =  l:interv 

index  =  [a  b  c  f  g  o]; 
muy=0; 


%(*  calculate  value  for  /x,  equation  (27)  *) 
for  il  =  l:nuser-l 

muy=muy+(r(l,index(l)/r(2,index(il  +  l)))/valpha; 
end 
muy  =  1  /sqrt(  1  +  muy) ; 

%(*  calculate  value  for  Pb,  equation  (25)  *) 
pb  =  (l-muy)/2 

%(*  calculate  the  value  for  PbM,  equation  (26)  *) 

pbm  =  0; 

for  cnt  =  0:m-l 

pbm =pbm  +  (comb((m-l  +  cnt),cnt)*(((l  +  muy)/2)/vcnt)); 
end 
pbm  =  (pb"m)*pbm; 

%(*  calculate  the  value  for  Pd  for  d  even,  equqation  (30)  *) 
for  q  =  l:6 

pd(q)  =  0.5*comb(d(q),d(q)/2)*pbm/Xd(q)/2); 
pd(q)  =  pd(q)*(l-pbmr(d(q)/2); 
for  e  =  l  +  d(q)/2:d(q) 

pd(q)=pd(q)  +  (comb(d(q),e)»(pbm/ve)*(l-pbm)^(d(q)-e)); 
end 
end 
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%(*  calculate  Pe,  equation  (29)  *) 
pe  =  ad*pd'; 
if  pe>  1,  pe  =  l;,  end 

%(*  calculate  the  integration  for  q(n),  equation  (38)  *) 
inc  =  (l-pe)"L*fr(l,o)*fr(2,a)*fr(2,b)*fr(2,c)*fr(2,f); 
inc=inc*fr(2,g); 
sum  =  sum+inc; 
end 
end 
end 
end 
end 
end 
sum  =  sum*delta(2)'\nuser-l)*delta(l); 
end 


Table  IX. 


SIMULATION    RESULTS    FOR    THE    PROBABILITY    OF 
CORRECT  PACKET  CAPTURE  FOR  RING  THREE 


q(l) 

q(2) 

q(3) 

q(4) 

q(5) 

q(6) 

M  =  l 

1 

0 

0 

0 

0 

0 

M=2 

1 

.0474 

0 

0 

0 

0 

M=4 

1 

.2578 

.0447 

.003 

.002 

0 

Plots  of  q(n)  versus  n  for  one-ring  and  three-ring  systems  are  shown  in  Figure 
11,  12,  13,  and  14. 
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Figure  11. 


Probability  Of  Correct  Packet  Capture  For  One-Ring  Network 
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Figure  12. 


Probability  Of  Correct  Packet  Capture  For  Ring  One 
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Probability  Of  Correct  Packet  Capture  For  Ring  Two 
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Figure  14.  Probability  Of  Correct  Packet  Capture  For  Ring  Three 
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APPENDIX  C 


%(*  initialization  *) 

nuser  =  15; 

p0=0.1; 

pr= 0.002; 

q  =  [0  1  0.9042  0.8099  0.7197  0.6350  0.5580  0.4901  0.4615  0.4285]; 

%(*  steady  state  value  of  q(n)  *) 
for  i  =  ll:nuser+4 

q(i)= 0.3852; 
end 

%(*  set  N  =  nuser,  calculate  7rnm,  it  is  called  pi(n,m)  *) 
N  =  nuser; 
forn  =  0:N 
nl  =  n+l; 
for  m  =  0:N 
ml  =  m+l; 
if  m<n-l,pi(nl,ml)  =  0; 

elseif  m=  =n-l 
sum  1  =  0; 
for  i  =  0:n 

il=i+l; 

suml  =  suml  +  comb(n,i)*pr^i*(l-pr)/v(n-i)*q(il); 
end 
pi(nl,ml)  =  (l-p0r(N-n)*suml; 

elseif  m>  =n 
suml  =  0; 

tempi  =  comb(N-n,N-m)*(l-p0r(N-m-l)*p0^(m-n); 
for  i  =  0:n 

il  =  i+l 

term  =  comb(n,i)*(l-pr)/v(n-i)*pr; 

term  =  term*((l-pO)*(l-q(il  +  m-n))  +  (N-m)*pO*q(il  +  m-n+l)/(m-n+l)); 

suml  =  suml  +  term; 
end 
pi(nl,ml)  =  tempi  *suml; 
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%(*  calculate  nn\  it  is  called  pis(n)  *) 
forn=l:N 

nl  =  n+l; 

suml  =  0; 

for  i=0:n-l 
il=i+l; 
suml=suml  +  (pis(il)*pi(il,nl-l)); 

end 

temp =pis(nl-l)-suml; 

pis(nl)  =  (l/pi(nl,nl-l))*temp; 

if  pis(nl)  <  0,pis(nl)  =  0;,end 
end 

%(*  calculate  7rn,  it  is  called  pin  *) 
den=sum(pis); 
pin = pis/den; 
%(*  calculate  fn  *) 
for  n=0:N 
nl  =  n+l; 
temp=0; 
total =0; 
for  k  =  0:N-n 
kl=k+l; 
sum  1=0; 
for  i  =  0:n 

suml  =  suml  +  (comb(n,i)*prAi*(l-pr)A(n-i)*q(i  +  kl); 
end 

temp  =  comb(N-n,N-n-k)*pOAk(l-pO)A(N-n-k); 
term=suml*temp; 
total  =  total  +  term; 
end 

f(nl)  =  total; 
end 

%(*  calculate  s  *) 
suml=0; 
fori  =  0:N 
il  =  i+l; 

suml  =  suml  +  (pin(il)*f(il)); 
s(il)  =  suml; 
end 


67 


%(*  calculate  the  Comb(n,m)  *) 

function  y  =  corab(n,m) 

y=xfact(n+  l)/(xfact(m+  l)*xfact(n-m  +  !)); 

%(*  calculate  the  factorial,  it  is  call  fact(n)  *) 

function  y=fact(x) 

ifx==0 

y=i; 

else 

y=i; 

for  i  =  l:x 
y=y*i; 
end 
end 

%(*  the  variable  xfact  *) 

x=0:52 

fori  =  1:53 

xfact(i)=fact(x(i)); 
end 
global  xfact 
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